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Abstract 

Within the Color Glass formalism, we construct the wavefunction of a high 
energy onium in the BFKL and large- A'^c approximations, and demonstrate the 
equivalence with the corresponding result in the Color Dipole picture. We propose 
a simple factorization formula for the elastic scattering between two non-saturated 
"color glasses" in the center-of-mass frame. This is valid up to energies which are 
high enough to allow for a study of the onset of unitarization via multiple pomeron 
exchanges. When applied to the high energy onium-onium scattering, this formula 
reduces to the Glauber-like scattering between two systems of dipoles, in complete 
agreement with the dipole picture. 
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1 Introduction 



Over the last decade, two different formalisms — the Color Dipole Picture (CDP) [1-12] 
and the Color Glass Condensate (CGC) [13-24] — have been developed to study high 
energy scattering in QCD. Both formalisms aim at a description of unitarization effects 
within perturbative QCD. But the specific technical realizations are quite different, and 
so are also the corresponding physical pictures. From the point of view of perturbative 
QCD, these are both leading logarithmic formalisms, in the sense that they resum the 
radiative corrections which are enhanced by powers of Ins, with s the total invariant 
energy squared. But unlike the BFKL equation 25 , to which they both reduce in the 
linear, or single scattering, approximation ("single pomeron exchange"), these formalisms 
allow also for non-linear effects, like multiple scattering, which are responsible for the 
unitarization of the scattering amplitude. Still, the way how the BFKL physics and the 
unitarization effects are encoded differs substantially from one formalism to the other. 
Thus, although expected, the equivalence between these two descriptions is by no means 
obvious (at least, not beyond the linear, BFKL, approximation). It is the main purpose 
of this paper to demonstrate this equivalence in the most explicit way, namely by using 
the Color Glass formalism to rederive the picture of high-energy onium-onium scattering 
that has been originally obtained within the Color Dipole formalism. 

A hint towards such an equivalence comes already from the fact that the equations 
for the non-linear evolution of the scattering amplitude are rather similar, although not 
exactly the same, in the two formalisms. For the CGC formalism, these are the Balitsky 
equations^ ^Hl; which form an infinite hierarchy : With increasing energy, the original 
projectile — say, a quark-antiquark pair in a colorless state, or "color dipole" — may 
radiate a gluon, so the scattering amplitude for the qq pair is naturally coupled in the 
evolution equation to the corresponding amplitude for the qqg system, and so forth. 
Within the CDP formalism, Kovchegov has managed to obtain a closed equation [TT] . 
but only after making the additional assumption that the color charges inside the target 
are uncorrelated. Kovchegov's equation is formally similar to the first equation in the 
hierarchy by Balitsky, and may be viewed as an approximation to the latter, but clearly 
it cannot be equivalent to it (since a closed equation contains less dynamical information 
than an infinite hierarchy). It is not a priori clear whether this lack of equivalence is 
intrinsic in the two formalisms (CDP and CGC), or merely related to the additional 
assumptions introduced by Kovchegov in his derivation. 

Our subsequent analysis will show that the two formalisms are in fact equivalent, 
at least, for the problem of interest here (onium-onium scattering at high energy) and 

^ These equations have been originally derived by Balitsky |19| . within a formalism using the operator- 
product expansion of Wilson line operators near the light-cone. Subsequently, Weigert has shown 
that this infinite hierarchy of equations can be compactly summarized into a single functional equation. 
Within the CGC formalism, the first equation in the hierarchy by Balitsky has been explicitly derived in 
Ref. '17', and the complete equivalence with the functional equation by Weigert (as far as the evolution 
of Wilson line operators is concerned) has been demonstrated in Ref. [221 ■ 
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within a wide (but limited) kinematical regime, which includes the BFKL regime and the 
onset of unitarization (see below for details). This implies that the results of the exact, 
numerical, investigations of the CGC theory, which are currently under way should 
match exactly the previous Monte-Carlo simulations of the CDP, due to Salam 0, but 
differ in their details |2Z] from the solution to the Kovchegov equation. 

At this stage, it is useful to briefly discuss the two formalisms, in order to emphasize 
their similitudes and differences. In both approaches, the main ingredient is the con- 
struction of the light-cone wavefunction of an energetic hadron in the leading logarithmic 
approximation. This means that, in the wavefunction, one keeps only the small-x gluons, 
which form a high density system and multiply — when further increasing the energy — 
via the basic Lipatov vertex. 

Within CDP [HI IH El E] (see also Ref . for a recent review and more references) , 
this construction is carried out in the large-iVc limit and in the linear (or BFKL) approxi- 
mation. The large- A^c limit allows one to treat a gluon like a qq pair in a color octet state. 
Then, the emission of a gluon from a primary "dipole" (a qq pair in a colorless state) is 
interpreted as the original dipole splitting into two: each new dipole is made of the quark 
(antiquark) component of the primary dipole and the antiquark (quark) component of 
the emitted gluon. By iterating this elementary process, one obtains a description of the 
evolved dipole (or "onium") wavefunction as a system of dipoles. Still because of the 
large-iVc limit, one can neglect the interference between emissions from different dipoles: 
the dipoles emit gluons independently, resulting in a tree of dipoles. The linear approxima- 
tion means that one neglects the interactions among the emitted dipoles, so that, e.g., the 
dipole number density evolves according to the linear BFKL equation^. This puts a high- 
energy limit on the applicability of the dipole picture: When the dipole density, which 
grows like N{Y) ~ e"^"^ (with F ~ Ins being the rapidity, and uq = (41n2)asAi'c/7r), 
becomes so large that alN{Y) ~ 1, non-linear effects like dipole recombination become 
important, and are expected to lead to saturation 1^ 1^ ITHl IT^ ITUI ITHj. 

Thus, gluon saturation is not included in the CDP wavefunction, and most likely it 
cannot be accomodated in this formalism (at least, not in a systematic way), since the 
dipole-dipole cross-section ~ is formally of higher order in the large-A^c counting. 
But, as we briefly recall now, the non-linear effects responsible for unitarization can be 
accomodated, namely, they can be naturally summed up in the scattering amplitude. This 
is what makes this formalism more suitable for the study of the high-energy scattering 
than the conventional BFKL approach. 

Specifically, consider onium-onium scattering in the center-of-mass frame (Yi = Y2 = 
Y/2) at an energy which is low enough for alN{Y/2) <^ 1, but high enough for alN'^{Y/2) ~ 
1. (Since N{Y/2) is a large number, these conditions leave a rather large window.) The 
first condition means that we can ignore saturation effects in the wavefunctions of any 
of the incoming onia. In the second condition, a'j.N'^{Y/2) is the probability that a pair 

^This is consistent with the large-A'c counting, since in the BFKL evolution each power of is 
multiplied by a factor of Nc. 
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of dipoles — one from each onium — scatters with each other in the 2-gluon exchange 
approximation. Together with the BFKL evolution of the individual wavefunctions, this 
gives the onium-onium scattering amplitude in the "single pomeron exchange" (here, 
BFKL pomeron) approximation. But when Y is so large that a'^N'^{Y/2) ~ 1, the scat- 
tering is so strong that "multiple pomeron exchanges" (i.e., the simultaneous scattering 
of several pairs of dipoles from the two onia) become equally important. A crucial sim- 
plification, however, is that this multiple scattering refers only to different dipoles: in the 
kinematical window of interest, the probability that a single dipole undergoes multiple 
scattering is still suppressed, since proportional to a'^N(Y/2). Because of this, the mul- 
tiple scattering series can be explicitly summed up, as shown in Ref. leading to a 
scattering amplitude of a generalized Glauber type which satisfies unitarity. 

Note that the judicious choice of the frame has been essential for the validity of the 
previous arguments: If, instead of the center-of-mass frame, we were to choose, say, the 
rest frame of the second onium [Yi = Y, Y2 = 0), then, at rapidities large enough for 
the unitarization effects to be important, the saturation effects in the first onium wave- 
function would be important as well (since N{Y) ~ N'^{Y/2)), and the GDP formalism 
would not be applicable any more. We see that the distinction between unitarization and 
saturation is frame dependent, and in an asymmetric frame the two phenomena cannot 
be disentangled from each other. To summarize, by working in the center-of-mass frame, 
GDP provides a simple description of the onset of unitarization, while avoiding the intri- 
cacies of the non-linear quantum evolution. But this has the drawback that the physics 
of saturation cannot be studied directly, but only indirectly, via its effects on the (boost 
invariant) scattering amplitude. 

On the other hand, the GGG formalism [H [H 113 El El dl (see also Ref. [21 
for recent reviews and more references) is precisely intended to provide a description 
of the non-linear effects in the hadron wavefunction, in particular, of saturation. In 
this approach, the only restriction in the construction of the wavefunction is the lead- 
ing logarithmic approximation, in the spirit of which the small-x gluons are treated as 
the products of radiation from fast moving "color charges" (the partons with higher val- 
ues of x), whose internal dynamics is "frozen" by Lorentz time dilation (thus forming a 
"color glass"). The hadron wavefunction at small-x is then fully specified by giving the 
probability law, or "weight function" , for the spatial distribution of these color charges. 

When further decreasing x, new quantum fluctuations become effectively frozen, and 
must be included in the color source. This can be done via a perturbative QGD calculation 
in which non-linear effects are taken into account via the coupling between the quantum 
fluctuations and the classical color field radiated by the sources constructed in previous 
steps. The result of this calculation is a functional renormalization group equation (RGE), 
sometimes referred to as the JIMWLK equation^, which governs the evolution of the 
weight function for the color sources with increasing Y = Inl/x. 

''This stands for Jalilian-Marian, lancu, McLerran, Weigert, Leonidov, and Kovner, which are the 
authors of Refs. ^1 El 120] in which this equation has been proposed and constructed. 
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So far, only approximate solutions to this equation have been constructed, which are 
valid in restricted kinematical domains^. In the limit where the color field are weak and 
the non-linear effects become negligible — this corresponds to not so high energies, where 
the gluon density is still low — , the RGE equation has been shown [121 E] to reproduce 
the BFKL equation for the gluon distribution In the opposite regime at very high 
energies, where the fields are strong, the non-linear effects in the RGE tame the rise of the 
gluon distribution with \ jx ("saturation"), and lead to the formation of a high-density 
gluonic state — the color glass condensate — characterized by a hard intrinsic scale, the 
saturation momentum Qs, and by large occupation numbers, of order l/a^, for all gluonic 
modes with momentum less than or equal to Qs |TSl |2S1 • 

One reason why these previous approximations are really crude, is that they have 
been merely concerned with two-point functions (like the gluon density), while in reality 
the functional RGE is equivalent to an infinite hierarchy of equations for the correlation 
functions. (This encompasses, in particular, the hierarchy by Balitsky P^HH^D]-) While in 
the strong field regime at saturation it seems to be extremely difficult to go beyond the 
mean field approximation of Refs. fSl 12^] in the weak field regime, on the other hand, 
one can rely on perturbation theory to simplify the RGE and study the coupled evolution 
of the various n-point functions. In this perturbative regime, one expects the RGE to 
reduce to the BFKL evolution, and, in particular, to the color dipole picture at large Nc- 

This was expected, but never proven. In this paper, we shall fill in this gap by 
showing that, under the assumptions alluded to before — weak fields and large-iVc — , 
the RGE which describes the evolution of the onium can be solved in the dipole basis, 
with a result which is indeed equivalent to that of GDP. In fact, we shall find exactly 
that representation of the onium wavefunction that has been used by Salam in numerical 
simulations of high-energy onium-onium scattering 

To emphasize that such an equivalence is not a priori obvious, let us mention here 
a few technical differences between the two approaches: Since built in terms of colorless 
dipoles, the GDP formalism is automatically free of infrared singularities, but displays 
ultraviolet divergences (the probability to radiate dipoles of arbitrarily small size is ar- 
bitrarily large) which cancel in between "real" and "virtual" contributions to physical 
observables. By contrast, in the GGG formalism there are no ultraviolet divergences, but 
since the corresponding degrees of freedom are now colorful, there are apparent infrared 
singularities, which cancel only in the calculation of gauge invariant quantities. By solving 
the RGE in the dipole basis, we shall reformulate the GGG formalism (in the low density, 
or BFKL, regime) in such a way that infrared finiteness becomes manifest. 

But when comparing the two approaches, what is most interesting is the way they 
describe the unitarization of high-energy scattering. As already explained, the GDP can 
do that only in a symmetric frame, like the center-of-mass frame, which for a given total 
energy minimizes the importance of the saturation effects, which are not under control. 

^The results thus obtained are consistent with analytic [TTIISIIEIESI and numerical |31''34"35''36[I37| 

studies of the Kovchegov equation, with investigations of the BFKL dynamics above Qs |28..32..38] . and 
also with previous studies of saturation [281 1151 1117) . 
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Besides, even in such a frame, CDP cannot be used at arbitrarily high energies, since 
saturation becomes eventually important. By contrast, the CGC formalism, which allows 
for non-linear effects in the hadron wavefunction, has no such a high energy limitation. 
But its applications to scattering are generally conditioned by the use of an asymmetric 
frame, like the infinite momentum frame of the target, in which the projectile has a simple 
structure (e.g., a qq pair), for which we know how to write down the scattering amplitude. 
This is why rigorous applications of this formalism have been so far restricted to physical 
situations which lend themselves naturally to such an asymmetric description, like deeply 
inelastic scattering, or proton-nucleus collisions (for references see ^^)- 

In principle, the fact of using different frames is not an obstacle against comparing 
the two approaches: Since the scattering amplitude is boost-invariant, this might be very 
well computed in an asymmetric frame within the CGC formalism, and the result then 
compared to that obtained by Salam IH] within CDP. But this would require an exact 
calculation using the fully non-linear CGC wavefunction, which is not yet available, and 
at best could be computed numerically [211 ■ To allow for an explicit, and more insightful, 
analytic comparison, we find it convenient to formulate the scattering problem in a sym- 
metric way also within the framework of CGC. Loosely speaking, we shall reformulate the 
CGC approach in such a way to mimic the strategy of CDP. 

More specifically, we shall show that the elastic scattering between two non-saturated 
color glasses can be represented as the eikonal coupling between the color charge in one 
glass and the light-cone Coulomb potential radiated by the color charge in the other 
glass. (The extension of this formula to the general case where one, or both, of the color 
glasses is saturated is complicated by the fact that we do not know how to write down 
the coupling between an arbitrary distribution of classical color charge and a strong non- 
Abelian field.) This formulation has the same limitations, and also the same advantages, 
as the CDP — it includes multiple scattering, but is inconsistent with non-linear effects 
in the hadrons wavefunctions — , so it can be used too to study the onset of unitarization 
in the center-of-mass frame. The comparison between the two formalisms becomes then 
straightforward, and their equivalence can be explicitly proven: When the incoming color 
glasses are (non-saturated) onia, we shall find that the CGC theory reproduces exactly 
the Glauber-like expression for the scattering amplitude originally derived by Mueller j3] , 
within the operator formulation of CDP. 

The plan of this paper is as follows: In Sect. 2, we shall construct the CGC description 
of an elementary dipole, and of the dipole-dipole scattering in the two gluon exchange 
approximation. The dipole wavefunction that we shall derive here will represent the initial 
condition for the quantum evolution with Y to be discussed in Sect. 3. Specifically, in 
Sect. 3, we shall show that, in the BFKL approximation and the large- A^^^ limit, the RGE 
for the quantum evolution of the color glass can be solved in the dipole basis (provided 
the initial condition is a dipole too). We shall thus recover the onium wavefunction of the 
color dipole approach. Finally, in Sect. 4, we shall propose a factorization formula for the 
scattering between two non-saturated color glasses, and show that this reproduces the 
CDP formula for onium-onium scattering in the center-of-mass frame. 
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2 An elementary dipole as a color glass 



In this section, we shall show that, for the purposes of high-energy scattering, an elemen- 
tary dipole (i.e., a quark-ant iquark pair in a color single state) can be described as a color 
glass (i.e., a random distribution of classical color sources with a specific weight function). 
In the spirit of the eikonal approximation, we shall assume the quark and the antiquark to 
be pointlike "particles" which propagate at (nearly) the speed of light, and whose trans- 
verse coordinates (with respect to the propagation axis) are frozen: x for the quark, and 
y for the antiquark. Since both the number of color charges and their spatial distribution 
are completely fixed, we anticipate that the only source of randomness — which makes 
the "glassy" description natural — refers to the color degrees of freedom: The fact that 
the dipole is a color singlet means that any average over the dipole wavefunction must 
include an average over color. In more conventional calculations, this averaging is per- 
formed by taking the trace over the color matrices which enter the coupling of the quark 
or the antiquark to an "external probe". (We emphasize that, in a scattering problem, 
this color averaging should be performed already at the amplitude level, and not only in 
the cross-section.) Alternatively, as we shall see, this averaging can be formulated as an 
integral over a set of random variables, with the interpretation of "classical color charges" . 
In what follows, we shall be mainly interested in the lowest-order scattering processes, 
which proceed via two gluon exchange; in that case, the random variables can be taken 
as Gaussian. 

To remain as simple as possible, we shall consider the scattering between two dipoles 
in the two-gluon exchange approximation, for which the result is well known. Here, we 
shall rephrase this standard result as the collision between two color glasses, in a form 
which is suitable for further generalizations, like the inclusion of quantum evolution and 
multiple scattering. 

To start with, we shall consider a more general process, whose description is well 
established in the CGC formalism, and which encompasses the dipole-dipole scattering 
as a special case, as we shall shortly see: This is the scattering between an elementary 
dipole and a color glass. The following discussion will also give us the opportunity to 
recall the basic ingredients of the CGC approach, and fix some notations. 

Quite generally, the ^-matrix element for a head-on dipole-hadron collision can be 
computed in the eikonal approximation as: 



where ^^(a;) and V{y) are Wilson lines describing the scattering of the quark, or the 
antiquark, off the color field in the hadron, and the color trace divided by Nc is the 
average over color alluded to before. Furthermore, the brackets in the right hand side 
indicate the average over the hadron wavefunction. The CGC theory provides an explicit 
realization for this average, namely (see, e.g., for more details) : 
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where a°'{x~,x) is the hght-cone Coulomb field radiated by the color sources in the 
hadron, W[a] is a positive-definite functional which specifies the probability to find a 
given field configuration (the "weight function"), and: 

V^[a] = Pexp (ig / dx~ a"" {x~ , x)tA . (2.3) 



with P denoting path-ordering in x~ . Note that we are using light-cone vector notations, 
= (t ± z)/ a/2, and our conventions are such that the hadron is a right mover (it 

propagates in the positive z, or positive x~^, direction), while the dipole is a left mover 

(negative z, or positive x~). 

The expression fl2.2|) is written in a specific gauge, namely the covariant gauge, in 

which a"' is the only non-zero component of the field in the hadron [A'^ = S^'^a"'), and is 

time-independent (i.e., independent of x~^). The relation between this field and the color 

charge density in the hadron p"'{x~,x) is simply given by the two-dimensional Poisson 

equation: 

-\/la''{x-,x) = p"(x",a;). (2.4) 

But the quantity computed in Eq. ()2.2|) is actually gauge-invariant, since so are both the 
weight function W[a] and the scattering operator built from Wilson lines |24j . 

Eq. (j2.2j) holds in any frame in which the dipole rapidity is not too large, so that one 
can neglect gluon radiation in the dipole wavefunction. If Y is the total rapidity gap, with 
y = 1/hadron + |2/dipoie|, the precisc couditiou reads (recall that we are working in a leading 
logarithmic approximation) : as||/dipoie| 1- On the other hand, there is no restriction 
on the rapidity |/hadron of the hadron. If as?/hadron ~ 1, this means that the effects of the 
quantum evolution must be included in the weight function, which to this purpose must 
be a function of ?/hadron ~ ^- The evolution of the weight function W[a] = WV[q;] with Y 
will be discussed in the next section. Here, we are only interested in the simple situation 
in which the hadronic target (i.e., the "glass" in Eq. ()2.2j) ) is itself an elementary dipole. 
This entails several simplifications: 

First, the color field of a dipole is weak. This means that the fiuctuating field a in 
Eq. (j2.2p has typically small amplitudes, ga <^ 1, so we can expand the Wilson lines in 
powers of ga. The lowest non-trivial contribution to 5* is of order^ g"^, and is obtained 
by expanding the Wilson lines to quadratic order in ga. (The terms of order g, which 
would be linear in a°'t°', vanish after averaging with the gauge-invariant weight function, 
or, alternatively, after taking the color trace.) We have: 



V^[a] ^ ^ + y dx-a''{x-,x)t'' (2.5) 

- Y J J dy~oi\x', x)a\y-, x) [e{x- - y-)et^ + e{y' - x-)th''] . 

^We refer here to the powers of g which are exphcit in the expansion of the Wilson hnes for fixed a 
(with ga <C 1). In the final result for S, additional factors of g may arise from the evaluation of the 
correlation functions of a as in Eq. H2.2|l . 
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In forming S, Eq. (j2.1|) . we see that the ordering of the color matrices in x becomes 
irrelevant in the present approximation, because of the symmetry of the color trace: 




ab 



We thus obtain: 



S{x,y) ^ 1 



(2.6) 



where: 




dx a"'{x ,x) 



(2.7) 



is the effective color field in the transverse plane, as obtained after integrating over the 
longitudinal profile of the hadron. We thus see that, in this approximation, the longitudi- 
nal structure of the color field becomes irrelevant as well (this is, of course, correlated with 
the fact that the ordering in x~ is unimportant). To evaluate the average in Eq. ()2.fi|l . it is 
therefore enough to consider the reduced weight function which specifies the distribution 
of the field a°'{x) in the transverse plane alone. 

Eq. ()2.6p holds whenever the color fields in the target are weak. In this regime, and 
to lowest order in perturbation theory, the scattering amplitude T = 1 — S" is determined 
solely by the two-point function of the color fields, which in turn can be traded in the 
(unintegrated) gluon distribution. In particular, in the case where the target is itself a 
dipole, Eq. ()2.6|) must reproduce the standard result for the dipole-dipole scattering in 
the two-gluon exchange approximation. This condition constrains the weight function 
W[a] which describes a dipole as a color glass, and can be used to actually construct this 
weight function (which, in the present approximations, is simply a Gaussian). 

Here, we shall perform this construction via a slight detour, which will allow us to 
introduce an intuitive representation for the color charge of a dipole. With this aim, we 
shall take advantage of the fact that the final expression for S must be symmetric in the 
two dipoles. So, our first goal will be to rewrite Eq. ()2.6p in such a way to make this 
symmetry manifest. Note that, to the order of interest, Eq. (j2.6j) is equivalent to: 



which is formally the same as Eq. (j2.H) but without the path-ordering. This rewriting 
has naturally introduced the following expression for the color charge density of the left- 
moving dipole: 



The hat on pi is to remind that this is a color matrix, as opposed to the c-number charge 
density p° in the CGC formalism. 




(2.8) 




(2.9) 
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Clearly, a similar color charge density should be associated also with the right-moving 
dipole, with the quark at Xq and the antiquark at yo : 

p"^(z) = ge[5^'\z - xo) - 5^'\z - y,)]. (2.10) 
This dipole generates the matrix- valued field (cf. Eq. ()2.4p ) : 

al{z) = ge[A{z-x^)-A{z-y^)\ = g{z\x^,y^), (2.11) 
where A(a; — y) is the two-dimensional Coulomb propagator: 

/cPk e'^^'^^^v') 1 1 

The infrared regulator /i is needed to write down the propagator, but is harmless in the 
present context, as it cancels out in the difference Q of the two propagators in Eq. (j2.1H) . 
This is a typical infrared cancellation permitted by the colorless nature of the dipole. 

These considerations suggest that, for dipole-dipole scattering, Eq. ()2.8p should be 
equivalent to the following, manifestly symmetric^, formula : 

S{xM^o,yo) ^ i^tri®trHe^/'^'^^£(^)"«(^). (2.13) 



As compared to Eq. ()2.8|1 . the field a has now the explicit expression ()2.1H1 . and the 
average over a reduces to the color trace {l/Nc)tTR. (In Eq. ()2.13|1 . the color traces act 
separately in the color spaces of the first dipole and the second dipoles, respectively. To 
make that clear, one could use different notations, say, t^, for the color matrices 

spanning these two spaces.) 

Let us check that Eq. ()2.13|) reproduces indeed the expected result: After expanding 
this expression to lowest non-trivial order, i.e., to quadratic order in each of the two color 
charge densities p2 and p^, and performing the color averages, one obtains: 

S{x, y\xo, t/o) ~ 1 - ^ (N! - 1) [V{x, y\xo, yo)?, (2.14) 



where 



V{x,y\xo,yo) = A{x - Xq) - A{x - yo) - A{y - Xo) + A{y - yo) 

1 {x - yo)"^ {y - xo^"^ 

— In — 

A-K [x - xo)^ [y - yo 



We recognize in Eqs. ()2.14|) - ()2.15|) the standard result for the dipole-dipole scattering 
amplitude in the two-gluon exchange approximation (see, e.g., |8j). 

Eq. fl2.13p makes explicit the fact that, for a dipole, the "average over the hadron 
wavefunction" reduces to an average over color, here implemented as a color trace (sepa- 
rately for each dipole). In what follows, we shall rephrase this in the color glass formalism. 

To make this symmetry even more obvious, note that ^ (P z pl{z)a%^{z) = J (fzV'al{z)V^a%{z). 
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That is, we shall return to c-number densities and fields, and replace the color traces with 
functional averages over a1 and a^, like in Eq. (|2.2|) . (It is now understood that, in a 
more symmetric notation, the field variable a"' in Eq. ()2.2j) should be renoted as a^.) 
We thus introduce c-number random color charge distributions of the form: 

p%iz) ^ Q%[5^'\z - xo) - S^'\z - yo)], (2.16) 

and similarly for pi, where the "classical color charges" Q'^ and Q1 are Gaussian random 
variables with the following correlation functions (s = L,R) : 

(q:)q = o, {q:qi,)q = 6'^'6ss'X, ^^Sr/ ^^-^^^ 

Note that corresponds to gt"-, and the two-point function in Eq. ()2.17p has been chosen 
to match {l/Nc)tT{t"'t^). It is easy to check that, to order g"^, Eq. ()2.13|) is equivalent to 

Six,y\xo,yo) = (^^^1 'i'-Pli-)'^%i-)^ ^ (2.18) 

with (cf. Eq. Q) : 

= Q%,[A{z - xo) - A{z - yo)] = Q'iiQ{z\xo,yo). (2.19) 

Then, we replace the Q-average in Eq. ()2.18|1 by functional averages like in Eq. ()2.2j) . 
Specifically, Eq. ()2.18|1 is the same as: 

Six,y\xo,yo) = J DM Wo[an] J DK] W^o[«l] e'^I '^'^^'-l(^)^'-ni^) , (2.20) 

with the following color glass weight function for an elementary dipole: 

Wo[a] = N j n^Q'^ ^^P{-^^} 5[<^\^)-Q"Q{z\xo,yo)\ (2.21) 

In Eq. ()2.2H) . Ng = N^ — l, A/" is a normalization factor, chosen such that / D[a] lVo[«] = 1; 
and the functional 5-function is understood with a discretization of the transverse plane: 

5[a\z)-Q'^g{z\x,,yo)\ = HH S{a^iz) ~ Q^giz\xo,yo)) 

z a=l 

Clearly, the weight function PFo[q^] depends also upon the dipole transverse coordinates 
Xq and yo, but this dependence is suppressed in its notation, for simplicity. 

The interaction piece exp{i JcPz p1{z)a'}^{z)} in the previous formulae, Eqs. (j2.20j) 
or ()2.18|) . may be recognized as the eikonal coupling between the color charge density 
in one system and the field created by the color charge of the other system. In the 
Abelian case (i.e., for electromagnetic dipoles), this would be the exact coupling at high 
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energy, including multiple collisions. But in the non-Abelian case, this is correct only to 
order g"^ (corresponding to a single scattering), as derived above. This is already clear 
from the fact that, in QCD, this interaction term is not gauge invariant. In general, to 
describe multiple (eikonal) scattering in QCD one has to use path-ordered exponentials, 
like in Eqs. ()2.1|) - ()2.3|) . to account for the non-commutativity of the color matrices in the 
interaction vertices. Still, in Sect. 4 below we shall argue that a simple factorized formula 
like Eq. (j2.2(J|) can be used also for multiple scattering, but only in a symmetric frame 
and within a limited range of energies. 

To conclude this section, let us derive an alternative expression for the dipole weight 
function, Eq. ()2.21|) . which will be more useful for what follows. Namely, since the color 
field Q'^Q of a dipole is weak, of order g (cf. Eqs. (j2.17|) and (|2.19|) ). it is possible to 
expand the (5-functional in Eq. (|2.2ip to quadratic order in Q'^Q without loss of accuracy. 
Then, one can explicitly perform the average over the classical color charges, and thus 
deduce the following formula, which is our final result in this section: 

Wo[a]=AfS^l + ^Jd^uJ d^v g{u\xo, yo)g{v\xo, y,) ^^.^^^^^.^^^ } ^N- (2-22) 

3 BFKL evolution: From color glass to color dipoles 

The main advantage of using a color glass description for an elementary dipole is that 
one can rely on the whole machinery of the CGC formalism El QHl 1^ to study the 
evolution of the dipole wavefunction with increasing energy (or rapidity Y). As mentioned 
in the Introduction, this amounts to solving a functional renormalization group equation 
(RGE) for the weight function WV[a] with the initial condition Wo[«] given by Eq. ()2.2ip 
or ()2.22|) . In principle, this equation can be used to compute VFy [a] up to arbitrarily high 
y, including in the non-linear regime at saturation. In practice, however, the complicated 
structure of the general RGE prevents us from obtaining explicit solutions, except un- 
der very crude approximations , or through numerical simulations (which present 
their own difficulties, though) j^Hl- So far, the solution to the RGE has not been fully 
investigated not even in the weak field regime at not so high energies, where saturation 
effects are unimportant, and BFKL physics should apply. In particular, this is the regime 
in which it makes sense to compare the predictions of the RGE with the Color Dipole 
picture of Refs. jSl IH El El • Although in this regime the RGE simplifies considerably JT] 
(see also below), it remains a non-linear equation, which is still difficult to solve. 

Throughout this paper, we shall restrict ourselves to this weak field regime, that we 
shall further simplify by taking the laige-N^ limit, in the spirit of the dipole picture. 
Under these assumptions, we shall be able to construct an explicit representation for the 
solution pyy[tt] in terms of a system of dipoles which undergoes BFKL evolution. In this 
representation, the equivalence with the color dipole picture of the onium wavefunction 
will become transparent. 
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Specifically, we shall find that, under the weak field and large- A^'c approximations, the 
RGE can be solved with the following Ansatz (compare to Eq. (j2.2H) ): 



where the notations are as follows: P/v(xi, yi, X2, 1/2', ■ ■ ■', xn, yN\y) is the probability den- 
sity to generate a system of dipoles with given transverse coordinates (namely, {Xi, y^) 
for the ith dipole) via the quantum evolution of an original dipole with coordinates {xq, yo) 
through a rapidity interval equal to Y. (Of course, P/v depends also upon the original 
coordinates {xq, yo), and so does but this dependence is kept implicit, to simplify 

the notation.) The dipoles are characterized also by the respective color charges (Q" for 
the quark and —Q^ for the antiquark, like in Eq. ()2.10j) ). but their distribution in color 
factorizes from that in the transverse space, and is separately a Gaussian for each dipole. 
The ^-functional enforces the total field in the system to be precisely the field generated 
by the A^ dipoles (in a given color configuration). Finally, there is a sum over config- 
urations, which includes the average over color (performed separately for each dipole), 
the integral J dTjsf over the tranverse coordinates of the A^-dipole system, and the sum 
over A^. The probabilities Pn{Y) are determined by solving a linear system of coupled 
evolution equations which follow from the functional RGE after inserting the Ansatz p.ip 
for the solution. As we shall see, these are precisely the equations used by Salam jS] in 
his numerical construction of the onium wavefunction based on the GDP. 

3.1 The RGE in the BFKL approximation 

In this subsection, after briefly recalling the general structure of the RGE, we shall derive 
its weak field, or BFKL, approximation, and discuss the simplifications which occur in 
this limit. 

The RGE is a functional Fokker-Planck equation, that is, a second-order functional 
differential equation of the diffusion type. It reads [d (see also Ref. (T^ for an early 
version of this equation, and Refs. [201 1^ for alternative derivations using the approach 
pioneered by Balitsky [TH]) : 



where J = J (fx, and the kernel ri°'^{x, y) is a positive-definite and non-linear functional 
of a, upon which it depends via Wilson lines: 




N=l 






(3.2) 




(3.3) 



12 



In this equation, 



lC{x, y.z)^-—^ f % = VlAix - z)V\/\{y - z\ (3.4) 

(ivr)^ [x — zY\z — yY 

and V^^ and V are Wilson lines in the adjoint representation, given by Eq. (j2.3p with 
t°- — i> T°. It is worth recalling here that the right hand side of Eq. (j3.2j) has been obtained 
after combining "real" and "virtual" contributions to the quantum evolution: The "real" 
contribution is represented by rj, while the "virtual" one is generated as the functional 
derivative of r] with respect to a. 

Note the label Y on the field argument ay of the functional derivatives in Eq. ()3.2j) : As 
explained in jT3 El] , this specifies the longitudinal coordinate x~ at which the functional 
derivatives are to be taken. However, this prescription becomes irrelevant in the weak field 
regime of interest, in which the (lowest order) dynamics is sensitive only to the projection 
of the field in the transverse plane, as defined in Eq. ()2.7p . Indeed, when ga 1, one 
can expand the Wilson lines in Eq. (|3.3p in perturbation theory, and obtain, e.g., 

(1 - V^KY^ ^ ^9{a''{x) - a%z)) (T^)^^, (3.5) 

which yields the lowest-order perturbative approximation to rj, of order g"^ : 

V^\x,y) ^ g'{T^T'')^,l ^}C{x,y,z) [a'^ix) - a%z)] [a\y) - a'iz)]. (3.6) 

As anticipated, Eqs. ()3.5p - ()3.6|) are insensitive to the longitudinal structure of the field, 
and so are also the observables computed in this approximation, like the dipole scattering 
amplitude ()2.6j) . This reinforces our conclusion in Sect. 2 that, in the weak field regime, 
it is enough to work with the reduced weight function which is a functional of a°'{x) alone. 
This is the functional that we shall denote as in what follows. Correspondingly, 

the argument of any functional derivative will be interpreted as a°'{x). 

The disappearance of the longitudinal coordinate from the problem is the first sim- 
plification specific to the weak field limit. The second simplification is that the kernel t] 
becomes just quadratic in a, as manifest on Eq. (j3.6p . Clearly, even with this kernel, the 
RGE ()3.2j) remains non-linear, but the non-linearity is now considerably simpler than 
with the general kernel ()3.3p . This is best appreciated by inspection of the evolution 
equations satisfied by the n-point functions (a(l)a;(2) ■ ■ ■a(?T,))y, which are obtained as 
follows: Start with the general definition of a correlation function in the CGC formalism: 

{0)y = J B[a] Wria] 0[a], (3.7) 

with 0[a] = q;(1)q;(2) ■ • -a^n), then take a derivative with respect to Y, use Eq. (13.2^ for 
dWy/dY, and integrate twice by parts in the functional integral, to finally obtain: 

f ' ,-ix,y)'0[a]\ , (3.8) 



BY \2 7_ 6a-{x) ' ' 6a\y) 



13 



If 1] is the general kernel (j3.3p . the r.h.s. of Eq. (j3.8|) involves n-point functions 
with arbitrary n, as generated by the expansion of the Wilson lines. Thus, the general 
RGE is equivalent to an intricate hierarchy of coupled evolution equations^, which are 
generally non-linear and must be solved simultaneously (since, e.g., the evolution of the 
2-point function is coupled to that of all the n-point functions with n >2). By contrast, 
with the quadratic kernel in Eq. (j3.6p . the evolution equation for correlation functions 
are linear, and do not mix n-point functions with different number of fields n. So far, 
only the equation satisfied by the 2-point function has been considered in the literature 
[TBI El and shown to be equivalent to the BFKL equation. In what follows, we 
shall analyze directly the functional RGE with kernel ()3.6p — thus encompassing all the 
n-point correlations — and use the large-iVc limit^ to construct an explicit solution. This 
is the solution anticipated in Eq. (|3.1|) . 

Before we conclude this subsection, there is one more issue which needs to be clarified: 
the convergence of the integral over z which enters the kernel rj, cf. Eqs. fl3.3|) or ()3.6|) . 
Clearly, there are no short-distance singularities: the poles in }C{x, y, z) aX z = x or 
z = y are compensated by the field-dependent factors in these equations, which vanish 
linearly at these points. But at large distances z = \z\ ^ x, y, we have /C ~ l/^:^, which is 
not enough to guarantee the absence of long-range (or "infrared" ) singularities. Whether 
such singularities appear or not, depends also upon the nature of the operator 0[a], and 
upon the behavior of the n-point functions of a at large transverse separations (i.e., upon 
the properties of the weight function). 

For a generic weight function, infrared finiteness has been verified so far (on specific 
examples) only for gauge-invariant operators, like the 5'-matrix element in Eqs. (|2.H) 
or ()2.6p . However, we shall see below that, for the onium weight function in Eq. ()3.ip . 
infrared finiteness is automatically ensured for any operator (whether gauge-invariant or 
not), because of the rapid decay of the color field of a dipole. 

3.2 Quantum evolution: the first step 

Before attacking the full RGE at arbitrary rapidity Y, let us study the very first step in the 
quantum evolution of a dipole. That is, start with an elementary dipole with transverse 
coordinates {xq, yo) at Yq = and study its evolution under a rapidity increment dF, 
with OgdY -C 1. Our aim is to show that this evolution can be viewed as the splitting 
of the original dipole into two new dipoles. This elementary example will also allow us 
to introduce in a simple setting some of the technical manipulations that will be useful 
later, in the general case. 

Since the weight function (j2.2H) of an elementary dipole is characterized by a single 

^This hierarchy becomes somehow simpler if the equations are written for the correlation functions of 
the Wilson lines (rather than a) . The resulting equations are those originally derived by Balitsky IT^ . 

^Note that the large-iVc approximation cannot be implemented at the level of the kernel rj alone, but 
requires some detailed information about the color structure of the weight function. 
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non-trivial correlation function, namely, the 2-point function: 

{a''{x)a''{y))^ = g^Cp G{x\xo, yo)G{y\xo, yo), (3.9) 

(we have also used Cp = Ng/2Nc), it suffices to compute the change in this quantity in 
the first step of the evolution, or, equivalently, its derivative at y = 0. This is obtained 
by letting Y in the evolution equation obtained by replacing 0[a] —>■ a"'{x)a"'{y) in 
Eq. (j3.8|) with kernel (j3.6p . Simple algebra yields^° : 

^{«)y = ^ /(2/C.,4<-<)(a^-<) 

— }CxxzO!y{a'^ — al) — ICyyza'^{ay — a'^)^ ^, (3.10) 

where = a"-{x), JCxyz = ^{x,y,z), and the factor of N^. has been obtained as 
{T'^T'^) ab^'^'^ = NcS""^. The first term within the brackets, proportional to JCxyz, repre- 
sents the "real gluon" contribution to the evolution, while the other terms make up the 
"virtual" contribution, and have been generated when commuting one of the functional 
derivatives in Eq. ()3.2|) through 7]. 

Consider first the convergence properties of the above integral over z. As expected 
from the general discussion, there is no singularity at short distances : the three terms 
within the integrand are separately ultraviolet finite. To study the large distance behavior 
{z ^ x,y), it is convenient to group separately the terms involving a^, and those without 
it. The latter combine to (cf. Eq. (jS3)) : 

{alal)^{2ICxyz - ICxxz - ICyyz] = i«)y^ (^-H) 

where the z-dependent prefactor decays like 1/z^ at large z, and the ensuing integral is 
convergent. Note that the potentially troublesome terms behaving like 1/z'^ have can- 
celled in the linear combination within the braces (i.e., in between "real" and "virtual" 
contributions). But for the terms involving one or two factors of a^, there is no such a 
cancellation. For instance: 

^-y.(«>y cx 4(«>y fo^^ z:$>x,y, (3.12) 

which leaves the place for a potential infrared problem. The only way to avoid this problem 
is that the 2-point function (a^a")^ decreases sufficiently fast at large separations \z — x\. 
There is no reason to expect such a property to hold in general. But it does hold for the 
case of interest here, i.e., for Y = and the 2-point function in Eq. ()3.9p . This is so 
because the color field of a dipole is rapidly decreasing with the separation from the 
center of the dipole : 

g{z\xo,yo) ~ 77--%- for ^ > ^0, yo, (3.13) 

ZTT 2^ 

"'^^Note incidentally that, if Eq. H3.1Q(I is used to deduce the evolution equation for the scattering 
amplitude T = 1 — 5 in Eq. 1)2. 6|l . one finds the BFKL equation, as it should. 
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where Vq = Xq — y^. Thus, when z is very large, the function in Eq. (|3.12p decays hke 
1/ ^ which is rapid enough for the convergence of the integral in Eq. (j3.1(jp . 
To conclude, the following equation is well defined: 



\ ^ VIY 



^ X 



g^Cp j ^K^xyz{Qx~ Qz){Qy - Qz) 

- ^xxz Qy{Qx - Qz) - ^yyz Gx{Gy " ^z) |; (3-14) 



(with the shorthand notation = Q{x\xQ,yQ)), so it is meaningful to perform an inte- 
gration by parts over z in its right hand side. With this aim, we shall use the second 
equality in Eq. ()3.4|1 which shows that, when Kxyz multiplies a function which vanishes 
a.t z = X and z = y (as in the integrand in Eq. ()3.14j) ). it can be replaced by: 

iCxyz -^l^l [^(^ - ^)^iy - ^)] • (3-15) 

It is tempting to use this replacement to perform the integration by parts. However, if 
one does so, one generates terms like V|^^ which are ill defined because of singularities 
ai z = Xq and z = y^. For instance: 



"^IGz^^g, = (VL6;(2|a;o,yo))' = (VLA(2-a;o)-V;A(2-t/o))' 
(27r)2 {xo-zyiyo-zy 



(3.16) 



has uncompensated poles at z = Xq and z = y^. 

To avoid this problem, we shall use the following trick, whose physical significance 
should become clear later: We first identically rewrite Eq. ()3.4|1 as 

^^y^ = ^VL(A(a;-z)-A(a;-a;o))VL(A(2/-z)- A(t/-a^o)) 

+ \ VL(A(a; -z)- A{x - y,))W\[A{y - z) - A{y - yo)) 

- ^ Vl{Aix - xo) - Aix - yo))Vl{Aiy - xo) - A{y - yo)) . (3.17) 

This is indeed the same as Eq. fl3.4j) since the new terms included within the braces are 
independent of z, and thus give zero when acted on by the derivatives. This equation too 
can be written as a total derivative when multiplying a function which vanishes at z = a:; 
and z = y : 

^xyz — ^ ^l^Q{x\xQ,z)g{y\xQ,z) + g{x\z,yQ)g{y\z,yQ) - ^(a;|a;o, 2/o)^(l/|a3o, t/o)}, 

(3.18) 

since Q{x\xq,z) = A{x — xq) — A{x — z), so that Vlg{x\xo, z) = — V^A(a3 — z) = 
S^'^\x — z), etc. But as compared to Eq. ()3.15p . the equation above has the advantage 
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that the function within the braces vanishes, by construction, at z = Xq and z = y^, so, 
after integration by parts, it will compensate the singularities at these points. 

Specifically, after inserting Eq. ()3.18|) into Eq. ()3.14|) and performing the integration 
by parts, all the terms but one — namely {V^QzY — cancel out, and we are left with: 



d 



Y 



Y=0 



27r 



^g'^Cp j M{xQ,yQ,z)l^-g{x\xQ,yo)g{y\xQ,yQ) 
+ g{x\xQ, z)g{y\xQ, z) + g{x\z,yQ)g{y\z,yQ)y (3.19) 



where the "dipole kernel" 

M{xo,yo,z) = 77^7 3.20) 

has been generated as in Eq. (j3.16p . Note that the surviving term (which gave rise to the 
dipole kernel) comes from the "real" piece in Eq. (|3.1(jp alone. But although the "virtual" 
piece has not given any explicit contribution to the final equation ()3.19|) . its presence in the 
original equation ()3.10|) was nevertheless essential to ensure infrared finiteness, and thus 
permit the manipulations leading to Eq. ()3.19|) . The latter is manifestly infrared finite 
because of the rapid decay of the dipole kernel at large distances: Mix^.y^.z^ 
when z > Xq, Dq. 

Although equivalent to Eq. ()3.14|) . as shown by the calculations above, Eq. ()3.19|) has 
the advantage to shed more direct light on the physical mechanism behind the quantum 
evolution of the dipole: According to Eq. ()3.9|) . the quantity g'^CFg{x\xQ, z)g{y\xQ, z) is 
recognized as the 2-point function (a°ay)o of the color field generated by a dipole with 
the quark at Xq and the antiquark at 2, and similarly for the other terms within the braces 
in Eq. (|3.19p . Furthermore, {g'^Nc/2Ti)M.{xQ^yQ, z) is the probability per unit rapidity 
for the emission of a soft gluon with transverse position z from a dipole with coordinates 
[xQ^yo). Equivalently, this is also the probability for the splitting of the original dipole 
{xq, yo) into two new dipoles with coordinates {xq, z) and (z, yo), respectively jS]. Thus, 
Eq. (j3.19j) has a transparent physical interpretation, which is illustrated in Fig. ^ The 
change in the field-field correlator is due to the splitting of the original dipole, which brings 
in new contributions to the field from the produced dipoles, but subtracts a corresponding 
contribution of the decaying dipole. 

This evolution can be translated into a change in the weight function, which, after 
one step in the quantum evolution, takes the following form (compare to Eq. (j2.2ip ) : 

WdyH = Pi(dr) / \{dQ^e-^5[a\x)-Q^g{x\xo,yo)\ (3.21) 

a=l 

+ I P2{z\dY) I \[Y[dQte-^5[a\x)-Qig{x\x,,z)-Q''^g{x\z 

*J Z J X 10 1 



2/0 j 



i=l,2 a=l 



17 



Figure 1: Graphical illustration of Eq. 



with the probabilities: 



Pi{dY) 



1-dY 



9 




M{xo,yo,z), 



P2{z\dY) 



dV 
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M{xo,yo,z). 



(3.22) 



Note that the integral over z in the formula for Pi has logarithmic singularities at z = xq 
and z = Hq, and thus must be computed with an ultraviolet cutoff (a minimal distance 
p : \z — xq\ > p and |z — yol > p)- Such singularities are to be expected in the wave 
function — since there is an infinite probability to emit an arbitrarily small dipole — , but 
they cancel out in physical quantities, as manifest, e.g., on the r.h.s. of Eq. (|3.19p . which 
is free of any (ultraviolet or infrared) problem, as explained before. Also, they cancel in 
the condition of probability conservation, which at this stage reads: 



Eq. ()3.21|) has indeed the structure anticipated in Eq. ()3.1|) . In the next subsection, 
we shall prove this structure for arbitrary Y (within the validity range of the weak field 
approximation) . 

3.3 The BFKL evolution of the onium weight function 

To discuss the general case, it is more convenient to use an alternative form of the weight 
function (j3.1|) . in which the average over color is explicitly performed. As we shall see, 
this new form is a generalization of Eq. ()2.22|) . but unlike the latter it cannot be obtained 
via the straightforward perturbative expansion of the 5-functional in Eq. p.ip . The 
reason is that, for our subsequent study of unitarization, we are interested in relatively 
high energies, where the average number of gluons in the onium wavefunction N{Y) is 
so large, N{Y) ~ l/cts, that it interferes with the perturbative expansion. Note that 
this is not in contradiction with the weak field assumption that we have used so far: For 
the strong field effects, like saturation, to be important, the dipole number should be 
even larger^^, namely, N{Y) ~ (Indeed, the probability that a given dipole within 

-'^-'^This condition can be also stated in terms of the strength of the typical fluctuations of the field a : 
the strong field regime corresponds to ga ^ 1 (cf. Sect. 3.2), or, for the 2-point function: as{aa) ~ f. 




(3.23) 
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the onium interact with any one of the other dipoles is of order a'^N{Y).) Still, already 
when asN{Y) ~ 1, the perturbative expansion needs to be reorganized, to permit the 
resummation of the terms of order {asN{Y))^ with n > 1. 

To illustrate the difficulty, consider just the second-order term in the expansion of the 
5-functional in Eq. ()3.H) . After performing the average over color, this term is formally 
of order (with the shorthand notation Qi{x) = Q[x\Xi,yi)) : 



In reality, however, this contribution is of order UgN , because of the terms in the sum, 
and thus is of order one for a typical configuration with ~ ^(^) ~ l/'^s- Similarly, 
there are higher terms in this expansion which are of order (agA^)", and thus contribute 
to leading order too. In order to isolate these terms, and perform the color average for 
them, we shall use the integral representation of the (5-functional : 

a--Y,Q■Q^ = /"D[e]e'^^^"(^)°'"(^)e-'S"-«?/=^«''(^)^>(^). (3.25) 
i=i -I 



The Gaussian integrations over the color charges are now easily performed, to give: 

N 



N 

a" 

i=l 



For each dipole i, the exponent \{j^^QiY is truly a small quantity, of order A = g'^/2Nc, 
so it is legitimate to preserve just the second order term in each exponential. After this 
expansion, the functional integral over the auxiliary variables can be also performed : 



a 

i=l 



i=l 



-j:q^s^ )_ - ^^+2 /^^(^)m^) r^^. (3.27) 



Thus, the expression for the weight function that we shall use in the remaining part 
of this section, and which is equivalent to Eq. 1)3.11) to the accuracy of interest, reads: 



5[a] . (3.28) 



We would like to show that this expression is a solution to the RGE fl3.2|) with kernel 
()3.6|) . and deduce in the process the evolution equations for the probabilities Pj\fiY). To 
that aim, we shall work with the evolution equation in the form ()3.8p . where 0[a] is an 

Since {aa)Y ^ c(sN{Y), as it should be clear from Eq. H3.1|l . this gives again the condition a'j.N{Y) ~ 1. 
By contrast, when asN(Y) ~ 1, the field amplitude is of order one, and ga is truly perturbative. 
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arbitrary operator; this equation is convenient since it involves a functional integral over 
a, which facilitates the work with the functional derivatives in Eq. (j3.28|) . 

After inserting the expression ()3.6|) for t], the r.h.s. of Eq. ()3.8|) can be decomposed 
into two pieces: a "real" piece, in which both functional derivatives act on 0[a], and a 
"virtual" piece, which is generated when one of these derivatives act on rj. As for the 
example of the 2-point function discussed in Sect. 3.2, the "virtual" piece is important to 
cancel infrared singularities within the "real" piece, but it gives no explicit contribution 
to the final equation that we shall establish. Therefore, in what follows we shall restrict 
ourselves to the "real" contribution, which reads: 



■■5 

,2 



dY 2ti 

N=l 



N=l ^ 

^ r \ r 

B[a]6[a]ll\l + - Q^HG^iv) 



2 



X 



/ IC{x, y, z) [al - al] [a^ - at] |^ , (3.29) 



after performing some integrations by parts in the functional integral. 

The next step is to compute the action of the functional derivatives coming from the 
weight function. We shall do this in the large- A'^c limit, to simplify the color algebra. Be- 
cause of the presence of 6[a] in the functional integral, it is clear that two of the functional 
derivatives must act on the explicit quadratic form in a. There are two possibilities: (i) 
both these derivatives refer to the same dipole, and (ii) they refer to different dipoles. 
However, the second possibility is suppressed at large (since, unlike the first type of 
contribution, it does not provide a factor of Nc), and will be neglected in what follows. 
We thus obtain: 

d{0)y ^ ^ g^N, 



dY 271 

N=l 

N 



X 



N=l 
j^^l Jx,y,z 



\ f ^ , , . P 1 5'0\a 



5aiSai I Sa'^5a'^ 



(3.30) 



«=o 



Next, we shall perform an integration by parts in the integral over z in the second 
line of Eq. ()3.30|) . Note the similitude between this integral and that involving the "real" 
piece in Eq. ()3.14j) . Like in that case, the function Qi{z) = A{z — Xi) — A{z — yi) has 
logarithmic singularities at z = Xi and z = y^, which would be amplified by a too "brutal" 
integration by parts, based on Eq. ()3.15j) . Once again, this difficulty can be avoided by 
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using the analog of Eq. (|3.18|) (here, separately for each dipole) : 
1 



xyz 



g{x\xi, z)g{y\xi, z) + g{x\z, yi)g{y\z, yi) - g{x\xi, yi)g{y\xi, yi 



(3.31) 



where the expression within the braces vanishes at z = Xi and z = y^. Still as in Sect. 
3.2, the integration by parts generates the "dipole kernel" (for the ith dipole) : 

\2 



Mix„yi,z) = {V^g,{z)y 



[Xi - yi 



(3.32) 



(27r)2 {x,-zy{y,-zy 

The three terms within the braces in Eq. ()3.3H1 are then naturally combined with the 
functional derivatives acting on 0[a], to finally yield: 



d{0) 



Y 



dY 
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f;/dr.p.(r)f:n{i + K/ 6 

7V=1 i=l j^i ^ ^\Ju 



6 



7VW)7 



M{xi,yi,z) 



g{x\x„yi 



(3.33) 
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0\a] 



a=0 



Since the operator 0[a] is arbitrary, the equation above is equivalent to an equation 
for dWy/dY, which is most suggestively written for Wy+dY = Wy + {dWy / dY)dY : 



oo „ 

Wy+aY[a] = J2 dT^ P^{Y) 

AT 1 

2t\t ^ C -1 ^ 

„■ — 1 J z _ ■ — 1 



(3 



+ dY 



27r 



N 



gj{x) 



5a% 



En 

4 = 1 JT^J 



g^ix) 



5a% 



2n 



X I M{x,,y,,z) l + ^(^1^*'^)^) ^ + ^(/ ^^^\^^y'- 



After comparing this to Eq. ()3.28|1 . the effects of the evolution become transparent: When 
increasing the rapidity in one step, a given A^-dipole configuration can either survive as 
it is, but with a smaller probability, or evolve by radiating one soft gluon, which in the 
large-iVc limit is equivalent to the splitting of one of the original dipoles into a pair of 
new dipoles. To demonstrate that this evolution is indeed consistent with Eq. ()3.28j) . we 
still have to show that Eq. ()3.34j) is of the form: 



oo„ ^ \ \ f C 

Wy^^y\a\ = J2J Pn{Y + dr) n [l + 2 ( ^'^^^ fc^J 



i=l 
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5[a] 



(3.35) 



and thus deduce the evolution equations for the probabihties. 

For this purpose, it is convenient to adopt an economical labelling of the dipoles, 
which takes into account the specificity of the quantum evolution: A newly produced pair 
of dipoles involves just one more transverse coordinate in addition to the two coordinates 
of the parent dipole (see e.g. Fig.[H). Thus, if one starts with the dipole {xq, Uq) at y = 
and follows its evolution with increasing Y , then one needs one more coordinate, say a^i, 
to describe the ensuing two dipole system, two more coordinates, say Xi and X2., for a 
three dipole system, and so on. 

Therefore, a A^-dipole configuration generated by the evolution up to rapidity Y 
can be labelled by — 1 transverse coordinates aji, 0^2, . . . iX^^i (physically, these are 
the positions of the emitted gluons), and thus can be characterized by a probability 
density Pm{xi,X2, . ■ ■ ,x^_i\Y). (As before, the dependence upon the original coordi- 
nates (a^o, !/o) is kept implicit.) With this labelling, the coordinates of the A^ dipoles 
are: (a;o,a3i), {xi,X2), ... {xn^i^Hq), and the measure for the phase-space integration 
is simply: 

cLVm = d^a;id^£C2 • • • d^a^AT-i . (3.36) 

In particular, the iih dipole in the configuration is that with coordinates (tCj-i, Xi). Thus, 
we have to replace retrospectively {xi,yi) — > {xi_i,Xi) in all the previous formulae. 

We are finally in a position to check that Eq. ()3.34|) can be indeed brought into the 
form ()3.35p . and deduce that, for this to be possible, the probabilities Pn must obey the 
following recurrence formula: 



PNix^,...,XN_i\Y + dY) 



1=1 



l-dY^y J Mix,_,,x,,z) 



Pn{xi, . . .,Xn-i\Y) 



+ dF — — ^ M{xi^i, Xi+i, Xi)PN-i{xi, ...,Xi,..., xn-i\Y) . (3.37) 



i=l 



In the second line, a hat on the argument Xi in Pn-i means that the coordinate Xi is 
actually missing from the respective configuration of A^ — 1 dipoles (that is, the ith dipole 
in that configuration has coordinates (a;j_i, ajj+i)). 

The r.h.s. of Eq. ()3.37p is recognized as the sum of two terms: a loss term and a gain 
term. While the loss term, which describes the emission of one gluon from the original 
A^-dipole configuration, can be easily read off Eq. (j3.34j) . the gain term, on the other 
hand, which describes the formation of the A^-dipole configuration via the splitting of one 
dipole in an original configuration with only A^ — 1 dipoles, is more subtle, and can be 
recognized only after a judicious change of variables in the last term in Eq. ()3.34j) . This 
is explained in the Appendix. 

Clearly, the recurrence formula (|3.37p can be also written as a set of coupled evolution 
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equation for the probabilities P^'- 



dPNjY) 
dY 



27r2 




Pn{Y) 



+ 



27r2 



2 Pn-i{x,\Y), 



(3.38) 



(s^i-l ~ Xij'^i^Xij^i — Xi) 



where we have omitted all the obvious arguments. These equations must be solved with 
initial conditions which follow from Eq. ()2.22p . namely P/v(l^ = 0) = ^ati. 

It is straightforward to check that the probability is correctly conserved by the evo- 
lution ()3.38|1 (see the Appendix) : 



In fact, because of the simple structure of the equations (j3.38p . in which P/v is coupled 
only to Pn-1, the probability conservation is even more stringent, in the sense that it 
holds already for configurations with neighbouring numbers of dipoles: 



In this equation, which will be proven in the Appendix, dP^/ /dY is the loss term in 
Eq. dSISHl), while dP^^^l^/dY is the gain term in the corresponding equation for Pn+i- 

The coupled equations ()3.38p can be solved iteratively for successively higher proba- 
bilities: For = 1, one has a closed equation for Pi{Y); once this is solved, its solution 
Pi{Y) is inserted in the equation with A^ = 2 to give a closed equation for P2{Y), etc. 
But in order for these equations to be well defined, they must be supplemented with 
an ultraviolet cutoff: Indeed, written as it stands, the integral over z in the loss term 
is afflicted with logarithmic singularities due to the poles of the integrand a,i z = Xi 
and z = Xi^i. These singularities reflect the fact that we cannot forbid the radiation of 
dipoles of arbitrarily small sizes. Rather, we shall require a minimal size p for the radiated 
dipoles, which plays the role of an ultraviolet cutoff, and upon which the probabilities 
Pn depend logarithmically. (For instance, Pi{Y;p) represents the probability that, after 
the evolution through the rapidity interval Y, the onium wavefunction contains no more 
dipoles of size larger than p other then the parent dipole.) But this cutoff dependence 
cancels out in the calculation of any physical quantity, which involves a sum over all 
probabilities. This is clear, for instance, from Eq. (j3.33p . which governs the evolution of 
an arbitrary observable: in the r.h.s. of that equation, the integral over the coordinate 
z of the newly emitted gluon is free of short-range (and also long-range) singularities, 
since, e.g., Q{x\xi,z) when z ^ Xi. A more speciflc example will be given in the 



N=l 




(3.39) 




(3.40) 
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Appendix, where we shall construct the evolution equation for the dipole number den- 
sity (which turns out to be the BFKL equation), and demonstrate in the process the 
cancellation of ultraviolet singularities between the loss and gain terms. 

To summarize, the weight function ()3.28|) (or ()3.1|) ). with the probabilities Pn{Y) 
evolving according to the linear system of equations ()3.38|) . is the solution to the RGE 
in the weak field, or BFKL, regime and in the large-A^^c approximation. This includes 
the regime in which the total number of dipoles (or gluons) is as large as N{Y) ~ l/a^, 
which corresponds to a gluon distribution dA^/dY of order one. But this fails to apply at 
energies so high that N{Y) ~ (i.e., the gluon distribution is of order l/os), where 
the saturation effects are important. 

These are precisely the assumptions used in Refs. Ej to construct the onium 
wavefunction in the color dipole picture. Thus, it makes sense to compare the two pictures, 
and after doing so, it turns out that they are equivalent indeed. This is most easily seen 
by comparing the previous results in this section with the version of GDP used by Salam 
in his Monte-Carlo simulations of the onium wavefunction: In Ref.jS], Salam has used a 
recurrence formula equivalent to Eq. ()3.37|) to numerically construct the onium. 

4 Onium— onium scattering at high energies 

We now have all the ingredients necessary to study the onium-onium scattering within the 
color glass formalism, up to energies which are high enough for the unitarization effects to 
play a role. This requires putting together the onium weight function Wy [ot] that we have 
constructed in the previous section, and the factorized formula for the elastic scattering 
between two color glasses that we have proposed in Sect. 2. 

4.1 Symmetric scattering of two color glasses 

More precisely, in Sect. 2 we have shown that the symmetric formula ()2.20|) can be used to 
describe the low-energy scattering between two elementary dipoles. Here, we shall argue 
that this formula can be extended to high energies provided both color glasses remain 
non-saturated (that is, they remain in the weak field regime, as characterized in Sect. 3). 
Specifically, the generalization of Eq. (j2.2(J|) to high energies reads: 

Sy = j DM Wy.yian] j DM W^.M e^/'^'^^'^^^)^'"-^ , (4.I) 

where the two weight functions are computed in the BFKL approximation (cf. Sect. 3), 
and the rapidities Y — y and y should be such that the weak field condition is satisfied 
for the two incoming systems. This implies a frame-dependent upper limit on the total 
rapidity Y up to which Eq. ()4.1|) can be used. Clearly, the optimal choice is the center- 
of-mass (CM) frame, y = Y — y = Y /2, since this allows the highest value for Y before 
saturation effects start to be important in any of the two color glasses. This highest value 



24 



is determined by the condition alN{Y/2) ~ 1, which, together with the BFKL estimate 
for the gluon number : N{Y) ~ e'^o^ (with uoq = {4:\n2)asNc/TT), provides the following 
upper limit for the validity of Eq. ()4.1|) in the CM frame: 

^ — In ^ . (4.2) 

That Eq. ()4.H) describes correctly the single (BFKL) pomeron exchange approxima- 
tion, should be rather obvious: After expanding the exponential there to second order, 
one generates the 2-point functions of the color fields in each color glass, and these are 
well known to obey the BFKL equation ^HjEI- (For the particular case of onium-onium 
scattering, this will be verified in the Appendix.) What is, however, more interesting 
is that Eq. ()4.1|) can be trusted also beyond the single-scattering approximation. That 
is, the multiple scattering series, as encoded in the higher order terms in the expansion 
of the exponential, is consistently described by this equation up to energies well above 
the onset of unit arizat ion. We shall not prove this in general, but simply show that, 
for onium-onium scattering, Eq. ()4.1|) reproduces the correct result, including multiple 
pomeron exchanges, as originally obtained in the color dipole picture jlj. 

But one can nevertheless understand why, in this formulation, there is no inconsistency 
between having a linear approximation for the wavefunctions and keeping non-linear terms 
in the collision: The point is that, even if the color fields are indeed weak in the individual 
wavefunctions, in such a way that a'^N{Y/2) <ti 1, the scattering between the two systems 
can still be strong, since enhanced by the number of gluons (or dipoles) in both systems, 
and thus of order {asN{Y/2))'^. This becomes of order one already for Y ~ Yc/2, showing 
that there exists an interesting range of intermediate rapidities, namely, 

YJ2 <Y < i;, (4.3) 

within which saturation effects are still negligible, but the unitarization effects set in. Of 
course, the existence of such an intermediate regime is specific to the CM frame. If we 
were to work in an asymmetric frame, like the rest frame of one of the two hadrons, then 
unitarization effects and saturation effects in the wavefunction of the energetic hadron 
would start to be important at the same energy, namely, for Y ~ Yc/2. 

One may wonder, what happens if one, or both, of the color glasses are saturated. 
What would be the generalization of Eq. ()4.1|) to that case ? (This would be interesting, 
e.g., to study collisions with Y > Yc in the CM frame.) As we know, in the CGC 
formalism there is no difficulty of principle of dealing with a saturated wavefunction. 
This is determined by the functional RGE, which is explicitly known in the non-linear 
regime (although rather tedious to solve there). The true difficulty, which prevents us 
from extending the symmetric formula 1)4.11) to the non-linear regime at saturation, is 
that we do not know how to couple a generic distribution of classical color charges to a 
strong classical color field. That is, we do not know how to generalize the eikonal coupling 
in Eq. ()4.H) to the strong field regime in the non-Abelian case. 
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In fact, it is not clear whether such a generahzation exists, even in principle. Recall 
from Sect. 2 that we do know how to couple a simple projectile, like a color dipole, to a 
color glass condensate. This involves Wilson lines which describe the multiple scattering 
of each elementary constituent in the projectile off the color field in the target. But 
this requires a detailed information about the color matrix structure of each constituent, 
which may be difficult, if not impossible, to encode in an average description of the color 
glass type. As we shall see in the following calculations, this difficulty does not appear 
for onium-onium scattering in the rapidity range ()4.3p since an elementary dipole from 
one onium undergoes only single scattering off the color field in the other onium. 

One may finally observe that using a symmetric formula like Eq. ()4.1|) at arbitrarily 
high energies it not really necessary. Once one accepts to work with saturated wave- 
functions, one can very well use an asymmetric frame, in which the factorization of the 
S'-matrix is better under control (recall, e.g., Eqs. ()2.1|) - ()2.2|) ). The advantage of the 
formula ()4.1|) is precisely to extent the use of non-saturated wavefunctions up to energies 
well above the threshold for the onset of unitarization. 



4.2 Application to onium— onium scattering 

Let us now use the onium weight function in Eq. p.ip to compute the S'-matrix element 
for onium-onium scattering in Eq. ()4.H) . Eq. ()3.1|) is more convenient for this purpose 
(as compared to the other form of the weight function, Eq. ()3.28|) ) since the functional 
integrals over and ai can be immediately performed, with the result: 

Sy = J2 dV^P^{Y/2) J2 / dT^,P^,iY/2) /exp i$^$^Q-gp(^b) ) _-(4.4) 

Af=l N'=l \ ^ i=l j=l ' I Q,Q 

In this equation, Q1 are color charges for the right-moving dipoles, Q"- are the corre- 
sponding quantities for the A^' left-moving dipoles, and (cf. Eq. ()2.15|) ) : 

T^iAi) = V{x,^i,x,\yj^i,yj) (4.5) 

is a shorthand notation for the interaction potential between the right-moving dipole 
with coordinates {xi^i,Xi) and the left-moving one at (t/j_i,7/j). The brackets refer to 
the average over color, to be performed separately for each of the N x N' dipoles (cf. 
Eq. ()3.1|) ). It is straightforward to compute this average explicitely for one set of color 
variables, say, those associated with the right-movers. This gives: 



SnxN' ^ ^e^S-«"S.-«P(^l^)^^ = ^expj " ^ E [ E '^I ^^(^b) (4-6) 



The remaining average over Q is not so easy to perform, so we shall evaluate it by using 
various approximations. To understand the nature of these approximations, it is useful to 
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consider first the lowest order term in the expansion of the exponential, which corresponds 
to the single-scattering approximation (as obvious after comparing with Eq. (|2.14j) ): 



^ N , N' N' V 
^one-scatt ^ 1 - - ( 0,^0:^ | | m) \ _ 

i=l \ j=l m=l I Q 

i=l j=l 

Since X'^Ng ^ g'^/A at large Nc, this contribution to the scattering amplitude T = 1 — S 
is of order a'^NN'. After averaging with the BFKL weight functions, as in Eq. ()4.4j) . this 
yields a contribution of order {asN{Y/2))'^ ~ a^e"^"^ ("the single pomeron exchange"), 
which becomes of order one for Y ~ Yc/2, cf. Eq. ()4.2p . Thus, for rapidities Y ^ Yc/2, 
the single-scattering approximation breaks down, and some of the higher order terms in 
the expansion of the exponential in Eq. ()4.6|1 must be taken into account. 

Specifically, we need to include correctly the terms of order {a'^NN')^ for any n, since 
they contribute to leading order, but we can neglect terms with fewer powers of or A^' 
(for a given power of a^), since, e.g., a'^N <^ 1 in the regime of interest, cf. Eq. ()4.3p . 
The strategy to isolate such terms is similar to that used in the construction of the weight 
function in Sect. 3.3. Namely, we first rewrite the exponential in Eq. ()4.6|) as a product 
of A^ exponentials, one for each dipole, and then we keep only the second order term in 
the expansion of each factor in this product: 

\ i=l ^ j=l m=l J ' Q 

This is justified since the exponent which has been expanded out is of order a^A^' (after 
averaging over color), and thus is genuinely small. Physically, this means that for a given 
dipole in the right-moving onium, it is sufficient to consider its single scattering with 
any of the dipoles in the left-moving onium. This is similar to the discussion leading 
to Eq. ()2.6|) for dipole-hadron scattering: whenever the hadron is characterized by weak 
fields (as is the case here for both onia), the scattering of a single dipole can be computed 
in the two-gluon exchange approximation. We thus see that, in the present context, the 
multiple scattering consists in the simultaneous scattering of several pairs of dipoles from 
the two onia, while each individual dipole undergoes, at most, single scattering. 

But even after the expansion leading to Eq. (|4.8p . the color average over Q remains 
difficult to perform, and will be evaluated here only in the large- A^c limit. This requires 
some clarifications, since the A'c-counting turns to be quite different in the construction 
of the wavefunction and in the scattering problem. Recall that, the dominant effects that 
have been resummed in the construction of the weight function in Sect. 3 were terms 
of order [asNfY)"-, with n > 1. In the present discussion of high-energy scattering, 
the dominant effects that we are about to include are powers of a^NN' ~ a^e'^°^. Thus, 
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although their content is very different, the effects of order asN^Y and e'^"^ are really 
treated on the same footing, namely, as effects of order one. The large- A/'c approximation 
is then obtained by neglecting terms which are suppressed by negative powers of Nc as 
compared to these leading order terms. As we show now, the calculation of the color 
average in Eq. ()4.8|) greatly simplifies in this limit. 

Consider, for definiteness, the term quadratic in A in the expansion of the product 
in Eq. (j4.8j) . (The term linear in A has been already evaluated in Eq. (j4.7p .) This term 
describes the simultaneous scattering of two pairs of dipoles, or "two pomeron exchange" : 

^ ^ l<i</<Af ^ j,m=l j',m'=l ' Q 

Note that there are A^(A^ — l)/2 x (A^')^ terms altogether, so the contribution (j4.9|) is of 
order {a'^NN'Y, as expected. The A^^c-counting becomes transparent after averaging over 
color (cf. Eq. (ITT7|l ) : 

(p'jQm QfQL'^^ = ^'^\^Ng6jmSj'm' + Ng[6jj'6mm' + Sjm'Smj')^ (4-10) 

The two contributions within the braces — of order Ng and Ng, respectively — are 
illustrated in Fig. |H The contribution in Fig.|21a, which is dominant at large Nc, involves 
the independent scattering between two pairs of dipoles: the pair and the pair 

{l,m). Each collision brings in a factor Ng = N^ — 1 from the sum over colors. In the 
process in Fig. |21b, the sum over colors gets closed only after mixing all the four dipoles 
(/ — >■ j — > i — i> m — i> /), thus giving rise to only one factor of Ng. This second process will 
be neglected in the large-A^^c approximation. 

The generalization to the higher order terms in the expansion of the product in 
Eq. (j4.8|) is immediate: the large-A^^c limit is tantamount to performing the average over 
color independently for each of the factors in this product. This yields: 

i-^E[^(^w]> (4-11) 

i=l ^ j=l ^ 

which is recognized as the S'-matrix element for a right-moving system of A^ dipoles which 
scatter independently off the color field created by a left-moving system of N' dipoles. 
Note, however, that such an independence holds only for the scattering between given 
configurations of dipoles. After averaging over all the configurations as in Eq. (|4.4p . the 
scatterings get correlated with each other, because of the correlations included in the 
weight functions. It is remarkable that the unitarity constraint SnxN' < 1 is satisfied 
already for fixed configurations (and not only in the average), as obvious on Eq. ()4.11|) . 

To the accuracy of interest, Eq. (j4.1H) can be rewritten in the manifestly symmetric 
form: 

Sn.n' = exp ' (4-12) 

i=i j=i ' 
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a) 



b) 



Figure 2: The two types of two pomeron exchanges according to Eq. ()4.10|) . The dipoles 
i and I are right-movers, while j and m are left-movers. The label a or 6 on a gluon line 
denotes the color exchanged through that line. 



which is formally like a Glauber approximation: The multiple scattering series is re- 
summed as the exponential of minus the amplitude ()4.7|) for a single scattering. But, 
once again, this exponentiation holds only configuration by configuration. After averaging 
over all such configurations, the resulting S'-matrix for onium-onium scattering: 

Sy = J2 d^NPN{Y/2) J2 / dT,,,P^iY/2) exp - ^J^H [^('\j)T '(4-13) 

N=l N'=l ^ i=l j=l ^ 

differs significantly from a simple exponential of the "single pomeron exchange". As 
discussed in Refs. jH IHl IHZI, this difference has dramatical consequences in the high- 
energy regime where the ^'-matrix is very small. Eq. ()4.13|) coincides, as anticipated, 
with the formula proposed in the framework of the color dipole picture in Ref. and 
used for numerical studies of unitarization in Refs. [SHEl- 
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A Appendix 



In this Appendix we collect some of the technical details that have been omitted in the 
discussion of the quantum evolution of the onium weight function in Sect. 3. Specifically, 
we shall verify that the recurrence formula ()3.37j) follows indeed from the functional evolu- 
tion in Eq. ()3.34|1 . then we shall check the probability conservation in Eqs. ()3.39j) - ()3.4()j) . 
and finally we shall introduce the dipole number density and deduce the correspond- 
ing evolution equation. This turns out to be the BFKL equation, as expected from the 
corresponding analysis within CDP |3|. 

It turns out that the key technical step behind all the subsequent manipulations is 
a special change of variables within the integral over the transverse coordinates of the 
dipoles, with measure (j3.36p . Let us illustrate this with the derivation of the gain term 
in the r.h.s. of Eq. ()3.37|) . (The loss term there can be trivially inferred from Eq. ()3.34|) .) 

Start with the gain term corresponding to a given value of in the r.h.s. of Eq. ()3.34j) . 
This means that we are studying the evolution from a configuration with A^ dipoles at 
rapidity F to a configuration with A^+1 dipoles at rapidity F + dF. This term is rewritten 
here for convenience: 

OisdY J d^xid'^X2 . . ■d'^XN-id'^z Pj\f{xi, . . . ,X]\r_i\Y) 
X ^ ( Xj] J [x,^i,z][z,Xi\ M{ Xi—\i Xi^ Z 

), (A.l) 

i=l ^ jj^i ^ 



in compact notations where ag = g'^Nc/27T and, e.g. 



\Xi-l,Z 



(A.2) 



(The functional derivatives act on 6[a], which is omitted here, for simplicity.) Note that 
we are using the dipole labelling introduced above Eq. ()3.36|) : thus, as compared to the 
original notations in Eq. ()3.34p . we have here replaced {xi,yi) — > {xi_i,Xi). 
We would like to show that Eq. ()A.1|) can be equivalently rewritten as 



N+l 



d^ccid^iCs . . . d'^XN dPj^l^{xi, xn\Y) JJ [xi^i, Xi] , (A.3) 

i=i 

where dPj^^^ is the gain term in Eq. ()3.37|) with A^ —>■ A^ + 1. With this aim, take 
a particular term in the sum over i in Eq. ()A.1|) (say, the ith term), and change the 
integration variables as follows: 

Xi-^yi, ... , Xi_i Ui I , z ^yi, Xi ^ yi+i , . . . , x^^i y^. (A.4) 

This implies M{xi^i,Xi,z) —> M{yi-i,yi+i,yi) and 

PNixi, . . .,xn~i\Y) -> Pwivi, ■ ■ .,yi-i,yi+i, . . ■,yN\Y) = Pjv(yi, . . . . . . ,t/jv|y). 



30 



After this operation, Eq. ()A.1|) is rewritten as: 

„ N+l N 

as dY / d^t/i . . . d'^yN JJ [vi-i, Vi] ^ M{yi-i, Vi+i, yi)PN{yi, ...,yi,..., yN\Y), 

1=1 i=l 

which is indeed of the form ()A.3|1 with (after renoting yi as cCj) : 

N 

dPjf|i(a3i, . . .,xn\Y) = asdY^M{xi^i,Xi+i,Xi)PN{xi, ...,Xi,.. .,xn\Y) . (A. 5) 

i=l 

As anticipated, this is the gain term in Eq. (j3.37|l with N N + 1. 

To verify the probabihty conservation in its most stringent form, namely (cf. Eq. ()3.40j) ) 

j dr^ dPli\Y) + j dr^v+i dPi^li(r) = o, (a.6) 

start with Eq. ()A.5|) for dP^^^{Y), and perform the change of variables in Eq. ()A.4|) in 
reverse order, to deduce: 

JdTN+idPl^l^{Y) = asdY j d^x^ . . .d^x^-i Pn{xi, . . . ,Xn-i\Y)Y^^ J M{Xi^i,Xi, z), 

which is recognized indeed as — J dT^ dP^\ according to Eq. ()3.37|) . 

Let us finally introduce observables which are obtained by averaging over the trans- 
verse coordinates of the dipoles, and show how to construct evolution equations for them. 
An example is the dipole number density, defined as follows: 

oo „ N 

nY{x,y\xo,yo) = ^ dT n Pn{xi, ... ,xn-i\Y) ^ 6'^'^\xi_i - x)6^'^\xi - y). {A.7) 

N=l i=l 

This is the density of dipoles with the quark located at x and the antiquark at y produced 
after the evolution of an initial dipole {xq, y^) through a rapidity interval equal to Y . More 
generally, an observable of this type, which is not sensitive to the dipole color charges and 
fields, but only to their transverse positions, is computed as follows: 

oo „ 

0{Y) = dTMPN{Xl,...,XN-l\Y)ON{x,,...,XN-l), (A.8) 

7V=1 

and satisfies an evolution equation which is obtained by taking a derivative w.r.t. Y in 
the above equation and using Eq. (j3.38|) for dP^ / dY : 

)(0 Ptnia)- 



= dTN-g^ON{Xi,...,XN-l) + J dTN+l-^ON+liXi,...,XN)j, 
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where we have used the fact that the gain term vanishes for = 1. Via the same 
manipulations as before (namely, the change of variables in Eq. flA.4|) . but in reverse 
order), one readily obtains: 



Qy = dTNPNixi,...,XN_i\Y)^ M{xi_i,x„z 

N=l i=l ^ 

X l^- On{xi, . . . ,xn^i) + On+i{xi,. . . ,Xi^i,z,Xi,. . . ,xn-i)^- (A.IO) 

If On is the sum of terms (one for each dipole), so like in Eq. ()A.7|) . it is clear that 
the only such terms which survive in the evolution equation ()A.10|) are those associated 
with the dipole which has split in the course of evolution. All the other terms cancel in 
between the loss and gain contributions. 

Specifically, for the dipole number density ()A.7j) . Eq. (jA.lOj) gives: 

duY °° /" ^ /■ 

-g^ = dTNPNixi,...,XN-i\Y)^ M{xi_i,Xi,z) (A.ll) 

N=l i=l "^^ 

X I - S{xi^i - x)6{xi -y) + 5{xi^i - x)5{z - y) + 5{z - x)5{xi - ?/)|. 



After simple manipulations, this is rewritten as a closed equation for ny 
d 



gynYix,y\xo,yo) = I \ - M{x,y, z)nYix,y\xo,yo) (A.12) 



+ M{x, z, y)nY{x, z\xo, t/o) + M{z, y, x)nY{z, y\xo, yo) 
This must be solved with the following initial condition (cf. Eq. ()A.7|) with P/v ^ni) '■ 
no{x,y\xo,yo) = 6^'^\xq - x)6^'^\yo - y). (A.13) 

It is easy to see that Eq. ()A.12|) is free of both infrared and ultraviolet singularities. It is 
probably less obvious, but nevertheless true, that this equation together with the initial 
condition ()A.13|1 is equivalent to the following equation, which puts the evolution into the 
original dipole (a^o, yo) '■ 

nY{x,y\xo,yo) = as f M{xo,yo, z)^^ - nY{x,y\xo,yo) 



dV 



+ nY{x,y\xo,z) + ny(a;, y|z, ^/o)}- (A.14) 



This is recognized as the BFKL equation, and is the standard way to describe the dipole 
number evolution in CDP |21 H] • The equivalence between Eqs. ()A.12|) and ()A.14|) is easy 
to demonstrate once one realizes that the solution to any of these equations (with the 
initial condition ()A.13|) ) has the following symmetry property, which can be checked on 
the explicit solution to Eq. ()A.14|) . as found e.g. in Refs. 0E| : 

{x ~ yynYix,y\xo,yo) = {xq - yoY nYixo,yo\x,y) . (A.15) 
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Eqs. ()A.12|) and ()A.14|) describe the same physical process — the BFKL evolution of a 
system of dipoles — , but they approach this evolution from different ends. In Eq. ()A.12|) . 
the evolution proceeds via dipole splitting at the highest rapidity end: When Y gets 
increased by dY, it is the measured dipole {x, y) which gets created or annihilated. By 
contrast, in Eq. ()A.14|) the increase in rapidity is used to "push backwards" the original 
dipole (a^O) Vo) by an amount dY . Then, this dipoles undergoes an additional splitting, 
which appears as being earlier with respect to the dipole configuration at rapidity Y . 
Thus, from this perspective, the system grows via splitting at its lowest rapidity end. 

Note finally that, by using the definition ()A.7|) for the dipole number density, the 
onium-onium scattering amplitude in the single pomeron exchange approximation (i.e., 
the first non-trivial term in the expansion of the exponential in Eq. (j4.13|) ) can be written 
in the familiar form [21 IH E] : 

ry''°™''™"(a^o, 2/o|a3i,yi) = j d^xd^xd^yd^y ny_y(x, x|a;o, ?/o) 

X T2-si--(a;,a;|?/,i,) n^{y,y\x^,y,), (A.16) 
where {xo,yo) and {xi,yi) are the parent dipoles in the two onia, and 

T'-^'-°\x,x\y,y) ^ ^ [V{x,x\y,y)Y 

is the dipole-dipole scattering amplitude in the 2-gluon exchange approximation, cf. 
Eqs. flO^ - (jTT3|l . Eq. (|Xl4ll implies that the amplitude (r06|l satisfies the BFKL 
equation, as expected. 
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